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Abstract.  The  goal  of  this  paper  is  to  demonstrate  how  the  internal  symmetry  of  the  N-body  celestial-mechanics  problem  can 
be  exploited  in  orbit  calculation. 

We  start  with  summarising  research  reported  in  (Efroimsky  2002,  2003;  Newman  &  Efroimsky  2003;  Efroimsky  &  Goldreich 
2003)  and  develop  its  application  to  planetary  equations  in  non-inertial  frames.  This  class  of  problems  is  treated  by  the  variation- 
of-constants  method.  As  explained  in  the  previous  publications,  whenever  a  standard  system  of  six  planetary  equations  (in  the 
Lagrange,  Delaunay,  or  other  form)  is  employed  for  N  objects,  the  trajectory  resides  on  a  9(Al-l)-dimensional  submanifold 
of  the  12(A-l)-dimensional  space  spanned  by  the  orbital  elements  and  their  time  derivatives.  The  freedom  in  choosing  this 
submanifold  reveals  an  internal  symmetry  inherent  in  the  description  of  the  trajectory  by  orbital  elements.  This  freedom  is 
analogous  to  the  gauge  invariance  of  electrodynamics.  In  traditional  derivations  of  the  planetary  equations  this  freedom  is 
removed  by  hand  through  the  introduction  of  the  Lagrange  constraint,  either  explicitly  (in  the  variation-of-constants  method) 
or  implicitly  (in  the  Hamilton-Jacobi  approach).  This  constraint  imposes  the  condition  (called  “osculation  condition”)  that  both 
the  instantaneous  position  and  velocity  be  fit  by  a  Keplerian  ellipse  (or  hyperbola),  i.e.,  that  the  instantaneous  Keplerian  ellipse 
(or  hyperbola)  be  tangential  to  the  trajectory.  Imposition  of  any  supplementary  constraint  different  from  that  of  Lagrange  (but 
compatible  with  the  equations  of  motion)  would  alter  the  mathematical  form  of  the  planetary  equations  without  affecting  the 
physical  trajectory. 

However,  for  coordinate-dependent  perturbations,  any  gauge  dilferent  from  that  of  Lagrange  makes  the  Delaunay  system  non- 
canonical.  Still,  it  turns  out  that  in  a  more  general  case  of  disturbances  dependent  also  upon  velocities,  there  exists  a  “generalised 
Lagrange  gauge”,  i.e.,  a  constraint  under  which  the  Delaunay  system  is  canonical  (and  the  orbital  elements  are  osculating  in  the 
phase  space).  This  gauge  reduces  to  the  regular  Lagrange  gauge  for  perturbations  that  are  velocity-independent. 

Finally,  we  provide  a  practical  example  illustrating  how  the  gauge  formalism  considerably  simplifies  the  calculation  of  satellite 
motion  about  an  oblate  precessing  planet. 

Key  words,  celestial  mechanics  -  reference  systems  -  solar  system:  general  -  methods:  N-body  simulations  - 
methods:  analytical  -  methods:  numerical 


1.  Introduction 

1.1.  Prefatory  notes 

On  the  6th  of  November  1766  young  geometer  Giuseppe 
Lodovico  Lagrangia,  invited  from  Turin  at  d’Alembert’s  rec¬ 
ommendation  by  King  Friedrich  the  Second,  succeeded  Euler 
as  the  Director  of  Mathematics  at  the  Berlin  Academy. 
Lagrange  held  the  position  for  20  years,  and  this  fruitful  period 
of  his  life  was  marked  by  an  avalanche  of  excellent  results,  and 
by  three  honourable  prizes  received  by  him  from  the  Academie 
des  Sciences  of  Paris.  All  three  prizes  (one  of  which  he  shared 
with  Euler)  were  awarded  to  Lagrange  for  his  contributions  to 
celestial  mechanics.  Among  these  contributions  was  a  method 
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introduced  by  Newton  and  Euler  and  furthered  by  Lagrange 
for  his  studies  of  planet-perturbed  cometary  orbits  and  only 
later  applied  to  planetary  motion  (Lagrange  1778,  1783,  1808, 
1809,  1810).  The  method  was  based  on  an  elegant  mathemati¬ 
cal  device,  the  variation  of  constants  emerging  in  solutions  of 
differential  equations.  This  approach  was  pioneered  by  Newton 
in  his  unpublished  Portsmouth  papers  and  described  very  suc- 
cintly  in  Cor.  3  ans  4  of  Prop.  17  in  Book  I  of  his  “Principia”. 
The  first  attempts  of  practical  implementation  of  this  tool  were 
presented  in  a  paper  on  Jupiter’s  and  Saturn’s  mutual  distur¬ 
bances,  submitted  by  Euler  to  a  competition  held  by  the  Lrench 
Academy  of  Sciences  (Euler  1748),  and  in  the  treatise  on  the  lu¬ 
nar  motion,  published  by  Euler  in  1753  in  St.  Petersburg  (Euler 
1753).  However,  it  was  Lagrange  who  revealed  the  full  power 
of  this  approach. 
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Fig.  1.  Two  coplanar  ellipses,  having  one  common  focus,  are  assumed 
to  rotate  about  this  focus,  always  remaining  within  their  plane.  Let 
a  planet  be  located  at  one  of  the  points  of  the  ellipses’  intersection, 
P,  and  assume  that  the  rotation  of  the  ellipses  is  such  that  the  planet 
is  always  at  the  instantaneous  point  of  their  intersection.  At  some  in¬ 
stant  of  time,  the  rotation  of  one  ellipse  will  be  faster  than  that  of  the 
other.  On  these  grounds  one  may  state  that  the  planet  is  rapidly  mov¬ 
ing  along  the  slower  rotating  ellipse.  On  the  other  hand,  though,  it  is 
also  true  that  the  planet  is  slowly  moving  along  the  faster  rotating  el¬ 
lipse.  Both  viewpoints  are  equally  valid,  because  one  can  divide,  in  an 
infinite  number  of  ways,  the  actual  motion  of  the  planet  into  a  motion 
along  some  ellipse  and  a  simultaneous  evolution  of  this  ellipse.  The 
Lagrange  constraint  (7)  singles  out,  of  all  the  sequences  of  evolving 
ellipses,  that  unique  ellipse  sequence  which  is  always  tangential  to  the 
physical  velocity  of  the  planet. 

Below  we  shall  demonstrate  that  the  equations  for  the  in¬ 
stantaneous  orbital  elements  possess  a  hidden  symmetry  sim¬ 
ilar  to  the  gauge  symmetry  of  electrodynamics.  Derivation  of 
the  Lagrange  system  involves  a  mathematical  operation  equiv¬ 
alent  to  the  choice  of  a  specific  gauge.  As  a  result,  trajec¬ 
tories  get  constrained  to  some  9-dimensional  submanifold  in 
the  12-dimensional  space  constituted  by  the  orbital  elements 
and  their  time  derivatives,  as  demonstrated  in  (Efroimsky  2002, 
2003;  Newman  &  Efroimsky  2003;  Efroimsky  &  Goldreich 
2003).  However,  the  choice  of  this  submanifold  is  essentially 
ambiguous,  and  this  ambiguity  gives  birth  to  an  internal  sym¬ 
metry.  The  symmetry  is  absent  in  the  two-body  case,  but  comes 
into  being  in  the  A-body  setting  (N  >  3)  where  each  body  fol¬ 
lows  an  ellipse  of  varying  shape  whose  time  evolution  contains 
an  inherent  ambiguity. 

Eor  a  simple  illustration  of  this  point  imagine  two  coplanar 
ellipses  sharing  one  focus  (Eig.  1).  Suppose  they  rotate  at  dif¬ 
ferent  rates  in  their  common  plane.  Let  a  planet  be  located  at 
one  of  the  intersection  points  of  these  ellipses.  The  values  of  the 
elliptic  elements  needed  to  describe  its  trajectory  would  depend 
upon  which  ellipse  was  chosen  to  parameterise  the  orbit.  Either 
set  would  be  equally  legitimate  and  would  faithfully  describe 
the  physical  trajectory.  Thus  we  see  that  there  exists  an  infinite 


number  of  ways  of  dividing  the  actual  planet’s  movement  into 
motion  along  its  orbit  and  the  simultaneous  evolution  of  the  or¬ 
bit.  Although  the  physical  trajectory  is  unique,  its  description 
(i.e.,  its  parametrisation  in  terms  of  Kepler’s  elements)  is  not. 
A  map  between  two  different  (though  physically  equivalent) 
sets  of  orbital  elements  is  a  symmetry  transformation  (a  gauge 
transformation,  in  physicists’  jargon). 

Lagrange  never  dwelled  upon  that  point.  However,  in  his 
treatment  (based  on  direct  application  of  the  method  of  varia¬ 
tion  of  constants)  he  passingly  introduced  a  convenient  math¬ 
ematical  condition  which  removed  the  said  ambiguity.  This 
condition  and  possible  alternatives  to  it  will  be  the  topics  of 
Sects.  1-3  of  this  paper. 

In  1834-1835  Hamilton  put  forward  his  theory  of  canoni¬ 
cal  transformations.  Several  years  later  this  approach  was  fur¬ 
thered  by  Jacobi  who  brought  Hamilton’s  technique  into  as¬ 
tronomy  and,  thereby,  worked  out  a  new  method  of  deriving 
the  planetary  equations  (Jacobi  1866),  a  method  that  was  soon 
accepted  as  standard.  Though  the  mathematical  content  of  the 
Hamilton- Jacobi  theory  is  impeccably  correct,  its  application 
to  astronomy  contains  a  long  overlooked  aspect  that  needs  at¬ 
tention.  This  aspect  is:  where  is  the  Lagrange  constraint  tacitly 
imposed,  and  what  happens  if  we  impose  a  different  constraint? 
This  issue  will  be  addressed  in  Sect.  4  of  our  paper. 

The  main  line  of  reasoning  and  the  principal  results  pre¬ 
sented  in  this  paper  are  the  following.  In  a  concise  introduction, 
to  be  presented  in  the  next  subsection,  we  derive  Eq.  (16)  which 
is  the  most  general  form  of  the  gauge-invariant  perturbation 
equation  of  celestial  mechanics,  written  in  terms  of  a  disturbing 
force.  Then  we  transform  it  into  (25),  which  is  the  most  gen¬ 
eral  gauge-invariant  perturbation  equation  expressed  through  a 
disturbance  of  the  Lagrangian.  The  next  step  is  to  show  that,  in 
the  case  of  velocity-dependent  perturbations,  the  equations  for 
the  Delaunay  elements  are,  generally,  not  canonical.  However, 
they  become  canonical  in  a  special  gauge  (which  is,  generally, 
different  from  the  customary  Lagrange  constraint). 

Our  other  goal  is  to  explore  how  the  freedom  of  gauge 
choice  is  supplemented  by  the  freedom  of  choice  of  a  coor¬ 
dinate  frame  in  which  to  implement  the  method  of  variation  of 
constants.  This  investigation,  carried  out  in  Sect.  3,  provides  a 
physical  example  of  a  nontrivial  gauge  being  instrumental  in 
simplifying  orbit  calculations.  In  this  way  we  demonstrate  that 
gauge  freedom  in  celestial  dynamics  is  not  only  a  cultural  ac¬ 
quisition,  but  can  be  exploited  to  improve  computations. 

1.2.  Osculating  elements  vs.  orbital  elements 

We  start  in  the  spirit  of  Lagrange.  Before  addressing  the 
A-particle  case,  Lagrange  referred  to  the  reduced  two-body 
problem, 

u  r 

r  +  4  -  =  0, 

J.Z.  y. 

r  =  r planet  “  r sun?  7^  =  6r(fUpianet  L  ^sun)?  (1) 

whose  generic  solution,  a  Keplerian  ellipse  or  a  hyperbola,  can 
be  expressed,  in  some  fixed  Cartesian  frame,  as 

r  =  /(Ci,...,C6,  0, 


r  =  g(Cu  ...,  Ce,  t) , 


(2) 
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where 


df 


C=const 


(3) 


Since  the  problem  (1)  constitutes  by  three  second-order  dif¬ 
ferential  equations,  its  general  solution  naturally  contains  six 
adjustable  constants  C,.  At  this  point  it  is  irrelevant  which  par¬ 
ticular  set  of  the  adjustable  parameters  is  employed.  (It  may 
be,  for  example,  a  set  of  Lagrange  or  Delaunay  orbital  ele¬ 
ments  or,  alternatively,  a  set  of  initial  values  of  the  coordinates 
and  velocities.) 

Following  Lagrange  (1808-1810),  we  take  /  as  an  ansatz 
for  a  solution  to  the  A-particle  problem,  the  disturbing  force 
acting  at  a  particle  being  denoted  by  AF': 


-  =  AF. 


r 


(4) 


Now  the  “constants”  become  time  dependent; 


r  =  /(Ci(f),-,C6(f),  f), 


(5) 


whence  the  velocity 

df  dt  ^  ^  dCi  df  ^""2^  dCi  df  ’ 

I  I 


(6) 


acquires  an  additional  term,  Yj(df  ldCi)(dCildt). 

Substitution  of  (5)  and  (6)  into  the  perturbed  equation  of 
motion  (4)  leads  to  three  independent  scalar  second-order  dif¬ 
ferential  equations  which  contain  one  independent  parameter, 
time,  and  six  functions  C,(0-  These  are  to  be  found  from  the 
three  scalar  equations  comprised  by  the  vector  Eq.  (4),  and 
this  makes  the  essence  of  the  variation-of-constants  method. 
However,  the  latter  task  cannot  be  accomplished  in  a  unique 
way  because  the  number  of  variables  exceeds,  by  three,  the 
number  of  equations.  Hence,  though  the  physical  trajectory 
(comprised  by  the  locus  of  points  in  the  Cartesian  frame  and 
by  the  values  of  velocity  at  each  of  these  points)  is  unique, 
its  parametrisation  through  the  orbital  elements  is  ambiguous. 
This  circumstance  was  appreciated  by  Lagrange,  who  amended 
the  system,  by  hand,  with  three  independent  conditions. 


V  ^  ^^0 

2-j  dCi  df 


(7) 


and  thus  made  it  solvable^.  His  choice  of  constraints  was  mo¬ 
tivated  by  both  physical  considerations  and  mathematical  ex¬ 
pediency.  Physically,  the  set  of  functions  (Ci(f), ...,  C6(f))  can 


*  The  treatment  offered  by  Lagrange,  as  well  as  its  gauge-invariant 
generalisation  presented  in  (Efroimsky  2002,  2003;  Newman  & 
Efroimsky  2003),  addressed  only  the  case  of  a  coordinate-dependent 
perturbation  AF(r).  The  treatment  presented  in  this  paper  covers  dis¬ 
turbing  forces  AF(r,  r,  t)  that  are  arbitrary  vector-valued  functions  of 
position,  velocity,  and  time. 

^  Insertion  of  (5)  into  (4)  makes  Eq.  (4)  a  vector  equation  of  the 
second  order,  with  respect  to  functions  C,(f).  This  is  equivalent  to 
three  scalar  equations  of  the  second  order.  However,  these  may  be 
interpreted  as  three  scalar  equations  of  the  first  order,  written  for  six 
variables  C,(f),  Ci{t),  provided  we  amend  them  with  six  trivial  equa¬ 
tions  Ci(t)  =  dC,(f)/df.  This  gives,  for  each  particle,  nine  scalar  first- 
order  equations  for  twelve  quantities  C,(f),  C,(f).  Three  constraints 
will  then  make  this  system  fully  defined. 


be  interpreted  as  parameters  of  an  instantaneous  ellipse  (in  a 
bound-orbit  case)  or  an  instantaneous  hyperbola  (in  a  flyby  sit¬ 
uation).  Lagrange  found  it  natural  to  make  the  instantaneous 
orbital  elements  C,  such  that,  at  each  moment  of  time,  this  el¬ 
lipse  (hyperbola)  coincides  with  the  unperturbed  (two-body) 
orbit  that  the  body  would  follow  if  the  disturbances  were  to 
cease  instantaneously.  In  other  words,  these  ellipses  (hyperbo¬ 
lae)  are  tangential  to  the  physical  trajectory.  This  is  the  well 
known  condition  of  osculation,  which  is  mathematically  im¬ 
plemented  by  (7).  The  appropriate  orbital  elements  are  called 
osculating  elements.  This  way,  Lagrange  restricted  his  use  of 
the  orbital  elements  to  elements  that  osculate  in  the  reference 
frame  wherein  ansatz  (5)  is  employed.  Lagrange  never  explored 
alternative  options;  he  simply  imposed  (7)  and  used  it  to  derive 
his  famous  system  of  equations  for  orbital  elements. 

Such  a  restriction,  though  physically  motivated,  is  com¬ 
pletely  arbitrary  from  the  mathematical  point  of  view.  While 
the  imposition  of  (7)  considerably  simplifies  the  subsequent 
calculations  it  in  no  way  influences  the  shape  of  the  physical 
trajectory  and  the  rate  of  motion  along  it.  As  emphasised  in 
Efroimsky  (2002,  2003)  and  Newman  &  Efroimsky  (2003),  a 
choice  of  any  other  supplementary  constraint 


z 


2L 

dCi 


df 


=  0(Ci,  ...,  Ce,  t). 


(8) 


<[>  being  an  arbitrary  function  of  time  and  parameters  C,, 
would  lead  to  the  same  physical  orbit  and  the  same  velocities^. 
Substitution  of  the  Lagrange  constraint  (7)  by  its  generalisa¬ 
tion  (8)  would  not  influence  the  motion  of  the  body,  but  would 
alter  its  mathematical  description  (i.e.,  would  entail  different 
solutions  for  the  orbital  elements).  Such  invariance  of  a  phys¬ 
ical  picture  under  a  change  of  parametrisation  is  called  gauge 
freedom  or  gauge  symmetry.  It  parallels  a  similar  phenomenon 
well  known  in  electrodynamics  and,  therefore,  has  similar  con¬ 
sequences.  On  the  one  hand,  a  “good”  choice  of  gauge  often 
simplifies  solution  of  the  equations  of  motion.  (In  Sect.  3  we 
provide  a  specific  application  to  motion  in  an  accelerated  co¬ 
ordinate  system.)  On  the  other  hand,  one  should  expect  a  dis¬ 
placement  in  the  gauge  function  O,  owing  to  the  accumulated 
error  in  the  constants. 

A  simple  case  of  the  same  motion  being  described  by 
two  different  families  of  instantaneous  ellipses  is  presented  in 
Eig.  1 .  Two  coplanar  ellipses  have  a  common  focus  and  are  ro¬ 
tating  about  it,  in  the  same  plane  but  at  different  rates.  Assume 
that  one  ellipse  is  rotating  rapidly  and  another  slowly.  Then  it 
will  be  legitimate  to  state  that  the  point  of  their  intersection,  P, 
is  moving  rapidly  along  a  slowly  rotating  ellipse.  At  the  same 
time,  it  will  be  right  to  say  that  it  is  moving  slowly  along  a 
swiftly  rotating  ellipse. 

Derivation  of  the  conventional  Lagrange  and  Delaunay 
planetary  equations  by  the  variation-of-constants  method  incor¬ 
porates  the  Lagrange  constraint  (Brouwer  &  Clemence  1961). 


^  In  principle,  one  may  also  impart  ®  with  dependence  upon  the 
parameters’  time  derivatives  of  all  orders.  This  would  yield  higher- 
than-first-order  time  derivatives  of  the  Ci  in  subsequent  developments 
requiring  additional  initial  conditions,  beyond  those  on  r  and  r,  to  be 
specified  in  order  to  close  the  system.  We  avoid  this  unnecessary  com¬ 
plication  by  restricting  O  to  be  a  function  of  time  and  the  C,. 
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Both  systems  of  equations  get  altered  under  a  different  gauge 
choice,  as  we  now  show.  The  essence  of  a  derivation  suitable 
for  a  general  choice  of  gauge  starts  with  (6)  from  which  the 
formula  for  the  acceleration  follows: 


d^r  dg  dg  dC,  dd> 

dfi  dt  ^  dCi  At  ^  At 

^  + V  ^  + 
d^t  -4-^  dCi  At  At 

I 


(9) 


Together  with  the  equation  of  motion  (4),  it  yields: 

^  +  iiL  +  y  + 

dfi  r  -4-^  dCi  df  At 

I 

r=\r\^  \f\. 


(10) 


The  vector  function  /  was  from  the  beginning  introduced  as 
a  Keplerian  solution  to  the  two-body  problem;  hence  it  must 
obey  the  unperturbed  Eq.  (1).  On  these  grounds  the  above  for¬ 
mula  simplifies  to: 


dg  dC,  dO 

44  dCi  At  At 

I 


(11) 


This  equation  describes  the  perturbed  motion  in  terms  of 
the  orbital  elements.  Together  with  constraint  (8)  it  consti¬ 
tutes  a  well-defined  system  which  may  be  solved  with  respect 
to  ACijAt.  An  easy  way  to  do  this  is  to  use  the  elegant  trick 
suggested  by  Lagrange:  to  multiply  the  equation  of  motion 
by  dfjdCn  and  to  multiply  the  constraint  by  -dgjdCn-  The 
former  operation  results  in 

<9/  fy  dg  ACj 
dC„  dCi  df 

\  j  •' 


df_ 

dCn 


.  ^  df  A^ 
AF - - - , 

dCn  At 


(12) 


while  the  latter  gives 


dg 

'dC„ 


V  ] 


y  d£^  dy 
4-/  dCi  At 


dg 

'dCn 


O. 


Having  summed  these  two  equalities,  we  arrive  at: 


dC; 


y  [c.  c,i  =  2L  AF  -  itL  ii:  _  i2L 

4lj  ^  df  ac..  ac...  At  ac... 


df  dfl)  dg 


At  ac, 


dC„  At  dCn 


(13) 


(14) 


[Cn  Cj\  Standing  for  the  unperturbed  (i.e.,  defined  as  in  the 
two-body  case)  Lagrange  brackets: 


[C.  C_,1  =  4^  -  4L  A.. 

■'  dCn  dCj  aCj  dCn 


(15) 


It  was  selected  above  that  O  is  a  function  of  time  and  of  the 
“constants”  C,,  but  not  of  their  time  derivatives.  Under  this  con¬ 
vention,  the  above  equation  may  be  shaped  into  a  more  conve¬ 
nient  form  by  splitting  the  full  derivative  of  into  a  partial  and 
a  convective  one,  and  by  moving  the  latter  to  the  left-hand  side: 


dCn  dC, 


df  df  50 

dCn  dC„  at 


dg 

dCn 


O 


(16) 


This  is  the  most  general  form  of  the  gauge-invariant  perturba¬ 
tion  equations  of  celestial  mechanics.  In  the  Lagrange  gauge, 
when  the  O  terms  are  absent,  we  can  obtain  an  immediate  so¬ 
lution  for  the  individual  dC,7df  by  exploiting  the  well  known 
expression  for  the  Poisson-bracket  matrix  which  is  inverse  to 
the  Lagrange-bracket  one  and  is  offered  in  the  literature  for  the 
two-body  problem.  (Be  mindful  that  our  brackets  (15)  are  de¬ 
fined  in  the  same  manner  as  in  the  two-body  case;  they  contain 
only  functions  /  and  g  that  are  defined  in  the  unperturbed,  two- 
body,  setting.)  In  an  arbitrary  gauge  the  presence  of  the  term 
proportional  to  50/5Cj  on  the  left-hand  side  of  (16)  compli¬ 
cates  the  solution  for  dC,7df,  but  only  to  the  extent  of  requir¬ 
ing  the  resolution  of  a  set  of  six  simultaneous  linear  algebraic 
equations. 

To  draw  to  a  close,  we  again  emphasise  that,  for  fixed  in¬ 
teractions  and  initial  conditions,  all  possible  (i.e.,  compatible 
and  sufficient)  choices  of  gauge  conditions  expressed  by  the 
vector  function  O  lead  to  a  physically  equivalent  picture.  In 
other  words,  the  real  trajectory  is  invariant  under  reparametri- 
sations  permitted  by  the  ambiguity  of  the  choice  of  gauge.  This 
invariance  has  the  following  meaning.  Suppose  the  equations 
of  motion  for  C\,  ...,  Ce,  with  some  gauge  condition  O  im¬ 
posed,  render  the  solution  Ci(f),  ...,  C6(f)-  The  same  equations 
of  motion,  with  another  gauge  O  enforced,  furnish  the  solu¬ 
tion  C,(0  that  has  a  different  functional  form.  Despite  this  dif¬ 
ference,  both  solutions,  C,(0  and  Cfi),  when  substituted  back 
in  (5),  yield  the  same  curve  r{t)  with  the  same  velocities  r{t). 
In  mathematics  this  situation  is  called  a  fiber  bundle,  and  it 
gives  birth  to  a  1-to-l  map  of  C,(0  onto  Cft),  which  is  merely 
a  reparametrisation.  In  physics  this  map  is  called  a  gauge  trans¬ 
formation.  The  entire  set  of  these  reparametrisations  constitute 
a  group  of  symmetry  known  as  a  gauge  group. 

Just  as  in  electrodynamics,  where  the  fields  E  and  B  stay  in¬ 
variant  under  gradient  transformations  of  the  four-potential 
so  the  invariance  of  the  orbit  implements  itself  through  the 
form-invariance  of  expression  (5)  under  the  aforementioned 
map.  The  vector  r  and  its  full  time  derivative  r,  play  the  role 
of  the  physical  fields  E  and  B,  while  the  Keplerian  coordi¬ 
nates  Cl, ...,  Ce  play  the  role  of  the  four-potential  A'^. 

No  matter  whether  the  role  of  constants  C,  is  played  by 
the  six  Kepler  variables  (e,  a,  Mq,  a>,  Q.,  i )  or  by  some  six 
combinations  thereof  (like,  say,  the  Delaunay  set  (27)  or  the 
Jacobi  set  or  the  Poincare  set),  these  constants  are  always  called 
“orbital  elements”.  Only  in  the  case  when  the  condition  (7) 
of  Lagrange  is  imposed  will  these  orbital  elements  be  called 
“osculating”.  This  is  called  an  osculating  solution,  because  the 
particles’  positions  and  velocities  are  instantaneously  tangent 
to,  i.e.,  touch,  the  Keplerian  ellipses  (or  hyperbolae).  It  is  easy 
to  see  that,  by  virtue  of  the  differentiation  chain  rule,  the  oscu¬ 
lation  condition  (7)  remains  form-invariant  under  a  transition 
from  the  Kepler  set  of  elements  to  any  other  set.  This  is  most 
natural,  because  the  osculation  condition  is  purely  geometric 
one  and  does  not  depend  upon  our  preferences  in  the  choice  of 
convenient  parameters. 

A  comprehensive  discussion  of  all  the  above-raised  issues 
can  be  found  in  Efroimsky  (2003).  The  interconnection  be¬ 
tween  the  internal  symmetry  and  multiple  time  scales,  both 
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in  celestial  mechanics  and  in  a  more  general  ODE  context,  is 
addressed  in  Newman  &  Efroimsky  (2003). 


2.  Planetary  equations  in  an  arbitrary  gauge 

2.1.  Lagrangian  and  Hamiltonian  perturbations 

We  can  proceed  further  by  restricting  the  class  of  perturbations 
we  consider  to  those  in  which  AF  is  derivable  from  a  perturbed 
Lagrangian.  Such  a  restricted  class  of  disturbances  is  still  broad 
enough  to  encompass  most  applications  of  celestial-mechanics 
perturbation  theory.  This  happens  largely  because  we  shall  con¬ 
sider  Lagrangian  perturbations  dependent  not  only  upon  coor¬ 
dinates  but  also  upon  velocities.  It  will  enable  us  to  describe 
in  a  perturbative  manner  not  only  physical  but  also  velocity- 
dependent,  inertial  forces  emerging  in  non-inertial  frames  of 
reference.  This,  in  its  turn,  will  provide  us  with  an  opportunity 
to  work  in  the  coordinate  system  precessing  with  the  planet  of 
interest. 

Let  the  unperturbed  dynamics  be  described  by  the  unit- 
mass  Lagrangian  Foir,  r,f)  -  r  ^/2  -  U{r,  f),  canonical  mo¬ 
mentum  p  -  r  and  Hamiltonian  'Hotr,  p,  t)  -  p  ^/2  +  U(r,  t). 
The  disturbed  motion  will  be  described  by  the  new,  perturbed, 
functions: 

r  ^ 

Fir,  r,  t)  -  JU)  +  A£  -  —  -  U(r,  f)  +  A£{r,  r,  t),  (17) 


p  -  r  + 


dA£ 

dr 


and 


^(r,  p,t)^pr-j:^^  +  U  +  A'H, 


(18) 


(19) 


AFi  being  introduced  as  a  variation  of  the  functional  form,  i.e., 
as  AFi  =  Fiir,  p,  f)  -  Fioir,  p,  f).  The  Euler-Lagrange  equa¬ 
tions  written  for  the  perturbed  Lagrangian  (17)  will  give: 

r  =  -^+AF,  (20) 

or 

where  the  new  term 


AF  = 


dA£ 

dr 


dt  {  dr  j 


(21) 


is  the  disturbing  force.  Expression  (21)  reveals  that  whenever 
the  Lagrangian  perturbation  is  velocity-independent,  it  plays 
the  role  of  the  disturbing  function.  Generally,  though,  the  dis¬ 
turbing  force  is  not  equal  to  the  gradient  of  A£,  but  has  an  extra 
term  generated  by  the  velocity  dependence. 

Examples  in  which  a  velocity  dependent  A£  has  been  used 
in  a  celestial  mechanics  setting  include:  the  treatment  of  inertial 
forces  in  a  coordinate  system  tied  to  the  spin  axis  of  a  precess¬ 
ing  planet  (Goldreich  1965)  and  the  velocity-dependent  cor¬ 
rections  to  Newton’s  law  of  gravity  in  the  relativistic  two-body 
problem  (Brumberg  1992). 

Our  next  goal  will  be  to  employ  the  above  formula  (21)  to 
translate  the  gauge-invariant  perturbation  Eq.  (16)  from  the  lan¬ 
guage  of  disturbing  forces  into  that  of  Lagrangian  disturbances. 


2.2.  Gauge-invariant  planetary  equations 


When  the  expression  (21)  for  the  most  general  force  emerging 
within  the  Lagrangian  formalism  is  substituted  into  the  gauge- 
invariant  perturbation  Eq.  (14),  it  yields: 


dC  j 

Z  [C.  Cj] 


df  dA£ 
dCn  dr 


dCn  df  \  dr  j 


(22) 


[Cn  Cj]  being  the  Lagrange  brackets  defined  through  (15). 
Since,  for  a  velocity-dependent  disturbance,  the  chain  rule 

d_^_d_^^  dA£  dr 
dCn  dr  dCn  ^  dr  dCn 

^  ^  ^  ^23) 

dr  dCn  dr  5C„ 
takes  place,  formula  (22)  may  be  reshaped  into: 


dC  j 

Z  [c.  Q]  ^  = 


dA£ 

dA£ 

50 

dCn 

dr 

\dCn  df 

dCn]  \ 

5AX\ 
dr  /■ 


(24) 


Next  we  group  terms  so  that  the  gauge  function  ®  everywhere 
appears  added  to  d(A£)fdr,  and  bring  the  term  proportional 
to  dCj/dt  to  the  left  hand  side  of  the  equation: 


Z  {  [C«  Cj]  + 


df  d  ldA£ 
dCn  dCj  \  dr 


-H  <1> 


d 

dCn 


AX.- 


dg 


df  d 


dA£  d 


dCn  dCn  dt  dr  dC„ 


®  -H 


df 


5AX\ 
dr  )' 


(25) 


These  modifications  help  us  to  recognise  the  special  nature  of 
the  gauge  O  =  -d(A£)  j  dr,  which  will  be  the  subject  of  dis¬ 
cussion  in  the  next  subsection. 

Contrast  (25)  with  (16):  while  (16)  expresses  the  variation- 
of-constants  method  in  the  most  general  form  it  can  have  in 
terms  of  the  disturbing  force  AF(r,  r,  f),  Eq.  (25)  renders  the 
most  general  form  in  the  language  of  a  Lagrangian  perturba¬ 
tion  A£(r,  r,  f). 

The  applicability  of  so  generalised  planetary  equations  in 
analytical  calculations  is  complicated  by  the  nontrivial  nature 
of  the  left-hand  sides  of  (16)  and  (25).  Nevertheless,  the  struc¬ 
ture  of  these  left-hand  sides  leaves  room  for  analytical  simpli¬ 
fication  in  particular  situations.  One  such  situation  is  when  the 
gauge  is  chosen  to  be 

+t/(G  0,  (26) 

dr 

t](r,  f)  being  an  arbitrary  vector  function  linear  in  r.  (It  may 
be,  for  example,  proportional  to  r  or  may  be  equal,  say,  to 
a  cross  product  of  r  by  some  time-dependent  vector.)  Under 
these  circumstances  the  left-hand  side  in  (25)  reduces  to  the 
Lagrange  brackets.  The  situation  becomes  especially  simple 
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when  5AX/5r  happens  to  be  linear  in  r,  in  which  case  we  may 
put  t](r,  t)  -  dt^mdr  and,  thus,  employ  the  trivial  Lagrange 
gauge  0  =  0  instead  of  the  generalised  Lagrange  gauge.  We 
shall  encounter  one  such  example  in  Sect.  3.4. 

As  already  stressed  above,  the  Lagrange  brackets  are 
gauge-invariant  because  functions  /  and  g  are  defined  within 
the  unperturbed,  two-body,  problem  (1-3)  that  lacks  gauge 
freedom.  For  this  reason,  one  may  exploit,  to  solve  (25),  the 
well-known  expression  for  the  inverse  of  matrix  [C,  C j].  Its 
elements  are  simplest  (and  are  either  zero  or  unity)  when  one 
chooses  as  the  “constants”  the  Delaunay  set  of  orbital  variables 
(Plummer  1918): 


C,  =  {L,  G,  H,  Mo,  w,  Q) 

L  =  ■\/jIa,  G  =  -yJ/j-a  (1  -  e^),  H  =  -yj/ia  (1  -  e^)  cos  i,  (27) 


where  yu  s  G(msun  +  »Jpianet)  and  (e,  a,  Mq,  oj,  O,  i)  are  the 
Keplerian  elements:  e  and  a  are  the  eccentricity  and  semima¬ 
jor  axis.  Mo  is  the  mean  anomaly  at  epoch,  and  the  Euler  an¬ 
gles  w,  D,  i  are  the  the  argument  of  pericentre,  the  longitude 
of  the  ascending  node,  and  the  inclination,  respectively. 

The  simple  forms  of  the  Lagrange  and  Poisson  brackets  in 
Delaunay  elements  is  the  proof  of  these  elements’  canonicity 
in  the  unperturbed,  two-body,  problem:  the  Delaunay  elements 
give  birth  to  three  canonical  pairs  (Q,,  P,)  corresponding  to  a 
vanishing  Hamiltonian:  (L,  -  Mq),  (G,  -  cS),  (H,  -  Q).  In  a 
perturbed  setting,  when  only  a  position-dependent  disturbing 
function  R(r,  f)  is  “turned  on”,  it  can  be  expressed  through  the 
Lagrangian  and  Hamiltonian  perturbations  in  a  simple  manner, 
R(r,  f)  -  AX.(r,  f)  -  -A'7/(r,  f),  as  can  be  seen  from  the 
formulae  presented  in  the  previous  subsection.  Under  these  cir¬ 
cumstances,  the  Delaunay  elements  remain  canonical,  provided 
the  Lagrange  gauge  is  imposed  (Brouwer  &  Clemence  1961). 
This  long  known  fact  can  also  be  derived  from  our  Eq.  (25): 
if  we  put  0  =  0  and  assume  AX  velocity-independent,  we 
arrive  to 


dC  j 

Z  [C.  Cjl  ^ 


5  AX 

dCn  ’ 


(28) 


where 


AX  =  AX(/(C,  t),  t)^R(f(C,  t),  t)  = 

-A'HifiC,  t),  t).  (29) 

This,  in  its  turn,  results  in  the  well  known  Lagrange  system  of 
planetary  equations,  provided  the  parameters  C,  are  chosen  as 
the  Kepler  elements.  In  case  they  are  chosen  as  the  Delaunay 
elements,  then  (28)  leads  to  the  standard  Delaunay  equations, 
i.e.,  to  a  symplectic  system  wherein  the  pairs  (L,  -  Mo), 

(G,  -  ai),  (H,  -  Q)  again  play  the  role  of  canonical  variables, 
but  the  Hamiltonian,  in  distinction  to  the  unperturbed  case,  no 
longer  vanishes,  instead  being  equal  to  A'H  =  -  AX. 

In  our  more  general  case,  the  perturbation  depends  also 
upon  velocities  (and,  therefore,  AX  is  no  longer  equal 
to  -ATY).  Beside  this,  the  gauge  O  is  set  arbitrary.  As  demon¬ 
strated  in  Efroimsky  (2003),  under  these  circumstances  the 
gauge-invariant  Delaunay-type  system  is  no  longer  canonical. 
However,  it  turns  out  that  this  system  regains  the  canonical 


form  in  one  special  gauge,  one  that  coincides  with  the  Lagrange 
gauge  when  the  perturbation  bears  no  velocity  dependence.  The 
issue  is  explained  at  length  in  our  previous  paper  (Efroimsky  & 
Goldreich  2003).  Here  we  offer  a  brief  synopsis  of  this  study. 


2.3.  Generalised  Lagrange  gauge  wherein 
the  Delaunay-type  system  is  canonical 


Equation  (22)  was  cast  in  the  form  of  (25)  not  only  to  demon¬ 
strate  the  special  nature  of  the  gauge 


<D 


dA£ 

dr 


(30) 


but  also  to  single  out  the  terms  in  the  square  brackets  on  the 
right-hand  side  of  (25):  together,  these  terms  give  exactly  the 
Hamiltonian  perturbation.  Thus,  we  come  to  the  conclusion 
that  in  the  special  gauge  (30)  our  Eq.  (25)  simplifies  to 


dG  j 

Z  C/l  77 


dA'H 
'  dC„  ' 


(31) 


As  emphasised  in  the  preceding  subsection,  the  gauge  in¬ 
variance  of  definition  (15)  enables  us  to  use  the  standard 
(Lagrange-gauge)  expressions  for  the  matrix  inverse  to  [C„  C j], 
to  get  the  planetary  equations  from  (31). 

Comparing  (31)  with  (28),  we  see  that  in  the  general  case 
of  an  arbitrary  A£Xr,  r,  f)  one  arrives  from  (31)  to  the  same 
equations  as  from  (28),  except  that  now  they  contain  -A7T  in¬ 
stead  of  AX.  When  the  orbit  is  parametrised  by  the  Delaunay 
variables,  those  equations  take  the  form: 


dL 

dA'H 

d(-Mo) 

dAH 

(32) 

dr 

di-MoV 

dr 

dL  ’ 

dG 

dA'H 

d(-a)) 

dAH 

(33) 

dr 

d{-u)  ’ 

dr 

dG  ’ 

dH 

dAH 

d(-Q) 

dAH 

(34) 

dr 

5(-Q)’ 

dr 

dH 

which  is  a  symplectic  system.  Eor  this  reason  we  name  this  spe¬ 
cial  gauge  the  “generalised  Lagrange  gauge”.  In  any  different 
gauge  O  the  equations  for  the  Delaunay  variables  would  con¬ 
tain  O-dependent  terms  and  would  not  be  symplectic.  (Those 
gauge-invariant  equations,  for  both  Lagrange  and  Delaunay 
elements  are  presented  in  Efroimsky  2003  and  Efroimsky  & 
Goldreich  2003.)  This  analysis  proves  the  following. 


Theorem:  though  the  gauge-invariant  equations  for 
Delaunay  elements  are,  generally,  not  canonical,  they  be¬ 
come  canonical  in  the  “generalised  Lagrange  gauge”. 


That  this  Theorem  is  not  merely  a  mathematical  coincidence, 
but  has  deep  reasons  beneath  it,  will  be  shown  in  Sect.  4 
where  the  subject  is  to  be  approached  from  the  Hamilton- Jacobi 
viewpoint. 

The  above  Theorem  gives  one  example  of  the  gauge  for¬ 
malism  being  of  use:  an  appropriate  choice  of  gauge  can  con¬ 
siderably  simplify  the  planetary  equations  (in  this  particular 
case,  it  makes  them  canonical). 
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According  to  (18),  the  momentum  can  be  written  as: 


p  -  r  + 


dA£ 

dr 


=  Sr  +  <I)  + 


dA£ 

di- 


OS) 


which,  in  the  generalised  Lagrange  gauge  (30),  simply 
reduces  to: 


p^g.  (36) 

Vector  g  was  introduced  back  in  (2-3)  to  denote  the  functional 
dependence  of  the  unperturbed  velocity  upon  the  time  and  the 
parameters  C,.  In  the  unperturbed,  two-body,  setting  this  ve¬ 
locity  is  equal  to  the  momentum  canonically  conjugate  to  the 
position  r  (this  is  obvious  from  (18),  for  zero  AX).  This  way, 
in  the  unperturbed  case  equality  (36)  is  fulfilled  trivially.  The 
fact  that  it  remains  valid  also  under  perturbation  means  that, 
in  the  said  gauge,  the  canonical  momentum  in  the  disturbed 
setting  is  the  same  function  of  time  and  “constants”  as  in  the 
unperturbed,  two-body,  case.  Thus,  we  have  shown,  follow¬ 
ing  Goldreich  (1965),  Brumberg  et  al.  (1971),  and  Ashby  & 
Allison  (1993)"^  that  the  instantaneous  Keplerian  ellipses  (hy¬ 
perbolae)  defined  in  gauge  (30)  osculate  the  trajectory  in  phase 
space. 

Not  surprisingly,  the  generalised  Lagrange  gauge  (30)  re¬ 
duces  to  O  =  0  in  the  simple  case  of  velocity-independent 
disturbances. 


3.  Gauge  freedom  and  freedom  of  frame  choice 

3.1.  Osculating  ellipses  described  in  different  frames 
of  reference 


The  essence  of  the  variation-of-constants  method  in  celestial 
mechanics  is  the  following.  A  generic  two-body-problem  solu¬ 
tion  expressed  by 

r^fiC,  t). 


^  L 
dtJc  P  /’ 

is  employed  as  an  ansatz  to  solve  the  disturbed  problem: 
r^fidt),  t), 

df  df  dCi 
dt  dCj  df 


(37) 

(38) 

(39) 

(40) 

(41) 


dg  dg  dC,  dO 
dt  dCi  df  df 

_P_l  ^  ^  ^ 

P  f  ^  dCi  dt  ^  dt' 


(42) 


Treatment  presented  in  these  publications  was  equivalent  to 
choosing  the  generalised  Lagrange  gauge.  Below,  in  Sect.  3.3,  we 
shall  consider  that  development  from  the  viewpoint  of  gauge-invariant 
theory.  The  papers  Goldreich  (1965),  Brumberg  et  al.  (1971),  and 
Ashby  &  Allison  (1993)  are  unique  examples  of  non-osculating  el¬ 
ements  being  employed.  One  more  such  example  appeared  in  1987 
when  Borderies  &  Longaretti  (1987)  put  forward  their  theory  of  geo¬ 
metric  elements,  to  be  used  in  the  the  planetary  ring  dynamics. 


As  evident  from  (41),  our  choice  of  a  particular  gauge  is  equiv¬ 
alent  to  decomposition  of  the  physical  motion  into  a  movement 
with  velocity  g  along  the  instantaneous  ellipse  (or  hyperbola, 
in  the  flyby  case),  and  a  movement  associated  with  the  ellipse’s 
(or  hyperbola’s)  deformation  that  goes  at  the  rate  O.  It  is  then 
tempting  to  state  that  a  choice  of  gauge  is  equivalent  to  a  choice 
of  an  instantaneous  comoving  reference  frame  whereby  we  de¬ 
scribe  the  motion.  Such  an  interpretation  is,  however,  incom¬ 
plete.  Beside  the  fact  that  we  decouple  the  physical  velocity 
in  a  certain  proportion  between  g  and  O,  it  also  matters  which 
physical  velocity  (i.e.,  velocity  relative  to  what  frame)  is  decou¬ 
pled  in  this  proportion.  In  other  words,  our  choice  of  the  gauge 
does  not  yet  exhaust  all  freedom:  we  can  still  choose  in  which 
frame  to  write  ansatz  (40).  We  may  write  it  in  inertial  axes  or  in 
some  accelerated  system.  For  example,  in  the  case  of  a  satellite 
orbiting  an  accelerated  and  precessing  planet  it  is  convenient  to 
write  the  ansatz  for  the  planet-related  position  vector. 

The  above  kinematic  formulae  (40)-(42)  do  not  yet  contain 
information  about  our  choice  of  the  reference  system  in  which 
we  implement  the  variation-of-constants  method.  This  infor¬ 
mation  shows  up  at  the  next  stage,  when  expression  (42)  is  in¬ 
serted  into  the  dynamical  equation  of  motion  f  =  -(prlr^)  + 
AF  to  yield: 


dg  dCi  dO  5 AX  d  /5AX\ 

dCi  dt  dt  dr  dt  \  dr  ) 


(43) 


Complete  information  about  the  reference  system  in  which  we 
put  the  method  to  work  (and,  therefore,  in  which  we  define 
the  orbital  elements  C,)  is  contained  in  the  expression  for  the 
perturbation  force  AF.  For  example,  if  the  operation  is  car¬ 
ried  out  in  an  inertial  coordinate  system,  AF  contains  physical 
forces  solely.  However,  if  we  wish  to  implement  the  variation- 
of-constants  approach  in  a  frame  moving  with  a  linear  acceler¬ 
ation  a,  then  AF  also  contains  the  inertial  force  -a.  In  case  this 
coordinate  system  rotates  relative  to  inertial  ones  at  a  rate  p, 
then  AF  also  includes  the  inertial  terms  -2pxr-fixr-px 
(p  X  r).  In  considering  the  motion  of  a  satellite  orbiting  an 
oblate  precessing  planet  it  is  most  reasonable,  though  not  oblig¬ 
atory,  to  apply  the  method  (i.e.,  to  define  the  time  derivative) 
in  axes  that  precess  with  the  planet.  However,  this  reasonable 
choice  of  coordinate  system  still  leaves  us  with  the  freedom  of 
gauge  nomination.  This  will  become  clear  in  the  example  con¬ 
sidered  below  in  Sects.  3.2-3. 5. 


3.2.  Relevant  example 


Gauge  freedom  of  the  perturbation  equations  of  celestial  me¬ 
chanics  finds  an  immediate  practical  implementation  in  the  de¬ 
scription  of  test  particle  motion  around  an  precessing  oblate 
planet  (Goldreich  1965).  It  is  trivial  to  extend  this  to  account 
for  acceleration  of  the  planet’s  centre  of  mass. 

Our  starting  point  is  the  equation  of  motion  in  the  inertial 
frame: 


(44) 


where  U  is  the  total  gravitational  potential  and  time  deriva¬ 
tives  in  the  inertial  axes  are  denoted  by  primes.  Suppose  that 
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the  planet’s  spin  axis  precesses  at  angular  rate  /t(f)  and  that  the 
acceleration  of  its  centre  of  mass  is  given  by  a{f).  Physically, 
this  acceleration  originates  both  due  to  the  circumsolar  motion 
and  due  to  the  gravitational  pull  from  the  other  planets  (Krivov 
1993).  It  is  a  very  natural  (and  very  common  in  celestial  and 
planetary  dynamics)  technique  to  switch  to  a  comoving  or/and 
corotating  reference  frame,  in  order  to  better  visualise  the  prop¬ 
erties  of  the  physical  system.  For  example,  in  oceanography 
and  atmospheric  science  almost  all  work  is  carried  out  in  a 
frame  corotating  with  the  Earth.  In  our  case,  our  preference  will 
be  to  use  not  a  corotating  but  rather  a  coprecessing  frame,  i.e., 
a  coordinate  system  attached  to  the  planet’s  centre  of  mass  and 
precessing  (but  not  spinning)  with  the  planet.  In  the  new  coor¬ 
dinate  frame  the  inertial  forces  modify  the  equation  of  motion 
so  that  it  assumes  the  form: 

dU 

f  =  — - 2uxr-uxr-ux(uxr)-a,  (45) 

or 

time  derivatives  in  the  accelerated  frame  being  denoted  by  dots. 

To  implement  the  variation-of-constants  approach  in  terms 
of  the  orbital  elements  defined  in  the  accelerated  frame,  we  note 
that  the  disturbing  force  on  the  right-hand  side  of  (45)  is  gen¬ 
erated  according  to  (21)  by: 

AX  (r,  r,  t)  -  R  +  r  ■  (fix  r)  +  ^(fi  x  r)  ■  (fi  x  r)  -  a  ■  r,  (46) 

where  we  denote  by  R(r,  f)  the  gravitational-potential  pertur¬ 
bation  (which  is  the  perturbation  of  the  overall  gravitational 
potential  U).  Since 


5  AX 

dr 


-  fix  r. 


(47) 


the  corresponding  Hamiltonian  perturbation  reads: 


A-K  = 

-[R+p-(fixr) 


-4m 


a  -  r 


■]-  -[R  +  (rx  p)  ■  fi-  a  ■  r],  (48) 


with  vector  J  —  r  x  p  being  the  satellite’s  orbital  angular  mo¬ 
mentum  in  the  inertial  frame. 

According  to  (35)  and  (47),  the  momentum  can  be  writ¬ 
ten  as: 


variation-of-constants  Eq.  (43),  AX  is  given  by  formula  (46) 
and  A'H  by  (50). 

We  now  choose  to  describe  the  motion  in  the  generalised 
Lagrange  gauge  (30),  so  the  expression  (<^  +  fixr)  on  the 
right-hand  side  of  (50)  vanishes  (as  follows  from  (47)),  and  the 
expression  for  ATT  in  terms  of  /  and  g  has  the  form: 

A'H^-[R(f,t)+fi-(fxg)-a-f].  (51) 


At  the  same  time,  the  generic  expression  for  the 
variation-of-constants  method  given  by  (25)  simplifies 
to  (31).  Insertion  of  (51)  therein  leads  us  to 

[C,C,]^  =  ^  [R(f,t)+p-(fxg)-a-f].  (52) 

df  aCr 

Interestingly,  this  equation  does  not  contain  fi  even  though  it  is 
valid  for  nonuniform  precession. 

As  explained  in  Sect.  2.3,  in  the  generalised  Lagrange 
gauge  the  vector  g  is  equal  to  the  canonical  momentum  p  = 
r  +  d  AX.  I  dr.  In  the  case  when  the  velocity  dependence  of  AX 
is  called  into  being  by  inertial  forces,  the  momentum  is,  ac¬ 
cording  to  (47), 


p  -  r  + 


dAX 

dr 


r  +  fix  r, 


(53) 


which  is  the  particle’s  velocity  relative  to  the  inertial  frame  co¬ 
moving  with  the  accelerated,  rotating  frame.  In  this  sense  we 
may  say  that  our  elements  are  defined  in  the  accelerated,  rotat¬ 
ing  frame,  but  osculate  in  the  comoving  inertial  one. 

In  the  Appendix  we  provide  explicit  expression  for  each 
of  the  partial  derivatives  of  fi  ■  J  that  appears  in  the  planetary 
Eq.  (52). 


3.4.  Elements  defined  in  the  accelerated,  rotating 
frame,  that  osculate  in  this  frame 

Here  we  not  only  define  the  elements  in  the  accelerated,  rotat¬ 
ing  frame,  but  we  also  make  them  osculate  in  this  system,  i.e., 
we  make  them  satisfy  ®  =  0.  In  this  gauge,  expression  (50) 
takes  the  following  form: 

A-W  =  -  [R(f,  t)+fi-(fxg)  +  (iixf)-(iixf)-a-f],  (54) 


p^g  +  ^+fixf. 


(49)  while  Eq.  (25),  after  some  algebra^,  looks  like  this: 


whence  the  Hamiltonian  perturbation  becomes 
A'H^-[R  +  (fxg)-fi  +  (<^+fixf)-(fixf)-a-f].  (50) 

3.3.  Elements  defined  in  an  accelerated,  rotating 

frame,  that  osculate  in  the  comoving  inertial  frame 

In  this  subsection  we  recall  a  calculation  carried  out  by 
Goldreich  (1965),  Brumberg  et  al.  (1971),  and  Ashby  & 
Allison  (1993)  and  demonstrate  that  it  may  be  interpreted  as 
an  example  of  nontrivial  gauge  fixing. 

Let  us  implement  the  variation-of-constants  method  in  a 
frame  that  is  accelerating  at  rate  a  and  rotating  at  angular 
rate  fi  relative  to  some  inertial  system  S.  This  means  that,  in  the 


[C„  Q]  ^ 
at 


dAXf 


+fi  ■ 


(55) 


When  substituting  (54)  into  (55),  it  is  convenient  to  rent  the 
expression  for  ATT  apart  and  to  group  the  term  (px  f)  ■  (fix  f) 

^  Due  to  (47),  the  second  term  on  the  left-hand  side  in  (25)  is  pro¬ 
portional  to  [d(p  X  f)/dC  j\Cj  =  px  Q>  and,  therefore,  vanishes.  The 
second  term  on  the  right-hand  side  simplifies  in  accordance  with  the 
simple  rule  A  ■  (B  x  C)  =  (C  x  A)  ■  B. 
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with  the  last  term  on  the  right-hand  side  of  (55): 
at  dCn 


+R- 


df  ,  dg 

^xgr-/x  — 

uCfi 

df 


dC„ 


d 


dC„ 


In  so  writing  (56)  we  have  deliberately  cast  it  into  a  form  that 
eases  comparison  with  (52). 

In  the  Appendix  we  set  up  an  apparatus  from  which  the  par¬ 
tial  derivatives  of  the  inertial  terms  with  respect  to  the  orbital 
elements  may  be  obtained.  We  also  show  that  some  of  these 
derivatives  vanish.  However,  a  complete  evaluation  of  the  iner¬ 
tial  input  to  the  planetary  equations  in  the  ordinary  Lagrange 
gauge  involves  a  long  and  tedious  calculation  which  we  do  not 
carry  out. 


3.5.  Comparison  of  the  two  gauges 

One  of  the  powers  of  gauge  freedom  lies  in  the  availability  of 
gauge  choices  that  simplify  the  planetary  equations,  as  we  can 
see  from  contrasting  (52)  with  (56).  While  the  latter  equation 
is  written  under  the  customary  Lagrange  constraint  (i.e.,  for 
elements  osculating  in  the  frame  where  they  are  defined),  the 
former  equation  is  written  under  a  nontrivial  constraint  called 
the  “generalised  Lagrange  gauge”.  The  simplicity  of  (52),  in 
contrast  with  (56),  is  evident. 

By  identifying  the  parameters  C,  with  the  Delaunay  vari¬ 
ables,  one  arrives  from  (52)  and  (56)  to  the  appropriate 
Delaunay-type  equations  (see  Appendix  I  to  Efroimsky  & 
Goldreich  2003).  The  Lagrange  equations  corresponding  to 
(52)  and  to  (56)  may  be  easily  derived  from  each  of  these  two 
equations  by  choosing  C,  as  the  Kepler  elements  and  using  the 
appropriate  Lagrange  brackets. 

Although  the  planetary  equations  are  much  simpler  in  the 
generalised  Lagrange  gauge  than  in  the  ordinary  Lagrange 
gauge,  some  of  these  differences  are  less  important  than  oth¬ 
ers.  In  many  physical  situations,  though  not  always,  the  and 
fi  terms  in  (56)  are  of  a  higher  order  of  smallness  compared  to 
those  linear  in  ji,  and  therefore  may  be  neglected,  at  least  for 
sufficiently  short  times®. 


3.6.  Further  applications 

In  the  above  example  of  a  satellite  orbiting  a  wobbling  planet, 
an  evident  simplification  of  the  planetary  equations  (both  in  the 
Lagrange  and  Delaunay  forms)  was  achieved  through  imposi¬ 
tion  of  the  generalised  Lagrange  gauge.  This  optimal  gauge 
differed  from  the  standard  Lagrange  constraint,  because  in 
the  said  example  the  Lagrangian  perturbation  depended  upon 
velocities. 

A  similar  situation  emerges  in  the  relativistic  two-body 
problem.  In  the  relativistic  dynamics,  even  the  two-body 

®  As  an  example  of  an  exception  to  this  rule,  we  mention  Venus 
whose  wobble  is  considerable.  This  means  that,  for  example,  the  fi 
term  cannot  be  neglected  in  computations  of  circumvenusian  orbits. 


problem  is  perturbed,  the  role  of  perturbation  being  played 
by  the  relativistic  correction  to  the  Newton  law  of  gravity. 
Interestingly,  this  correction  depends  not  only  upon  the  posi¬ 
tions,  but  also  upon  the  velocities  of  the  binary  components 
(Brumberg  1992).  For  this  reason,  to  simplify  the  orbit  integra¬ 
tion  of  a  binary,  the  generalised  Lagrange  gauge  (not  the  cus¬ 
tomary  Lagrange  constraint)  should  be  imposed.  In  this  gauge 
the  calculations  will  be  very  considerably  simplified  (and,  for 
example,  it  is  in  the  generalised  Lagrange  gauge  that  the  equa¬ 
tions  for  the  Delaunay  elements  will  retain  their  canonicity). 

Another  simple  example  is  a  non-relativistic  reduced  two- 
body  problem  with  a  variable  mass.  In  this  case,  too,  the 
Lagrangian  acquires  a  velocity-dependent  correction.  Hence, 
in  this  case,  the  orbital  elements  will  be  convenient  to  intro¬ 
duce  in  the  generalised  Lagrange  gauge,  not  in  the  customary 
Lagrange  gauge. 


4.  Planetary  equations  and  gauges 
in  the  hamilton-Jacobi  approach 

In  this  section  we  demonstrate  that  the  derivation  of  planetary 
equations  in  the  A-particle  (N  >  3)  case,  performed  through 
the  medium  of  Hamilton-Jacobi  method,  implicitly  contains  a 
gauge-fixing  condition  not  visible  to  the  naked  eye.  We  present 
a  compact  account  of  our  study;  a  comprehensive  descrip¬ 
tion  containing  technical  details  may  be  found  in  Efroimsky 
&  Goldreich  (2003). 

The  Hamilton-Jacobi  analysis  rests  on  the  availabil¬ 
ity  of  different  canonical  descriptions  of  the  same  physi¬ 
cal  process.  Any  two  such  descriptions,  (q,  p,  'H(q,py)  and 
(Q,  P,  correspond  to  different  parametrisations  of 

the  same  phase  flow,  and  both  obey  Hamilton’s  equations.  Due 
to  the  latter  circumstance  the  infinitesimally  small  variations 

d6  =  pdq-^dt  (56) 

and 


dd^PdQ  -  Pi^dt 


(57) 


are  perfect  differentials,  and  so  is  their  difference 

-dW  =  d6-de^PdQ-pdq-('H*-P{)dt.  (58) 


Here,  vectors  q,  p,  Q,  and  P  each  contain  N  components. 
Given  a  phase  flow  parametrised  by  a  set  {q,  p,  P{(q,p,t)), 
it  is  always  useful  to  simplify  the  description  by  a  canonical 
transformation  to  a  new  set  (Q,  P,  p{*(Q,P,ty),  with  the  new 
Hamiltonian  Pt*  being  constant  in  time.  Most  advantageous  are 
transformations  that  nullify  the  new  Hamiltonian  P(*,  because 
then  the  new  canonical  equations  render  the  variables  (Q,  P) 
constant.  A  powerful  method  of  generating  such  transforma¬ 
tions  stems  from  (58)  being  a  perfect  differential.  It  is  sufficient 
to  consider  W  to  be  a  function  of  the  time  and  only  two  other 
canonical  variables,  for  example  q  and  Q.  Then  (58)  may  be 
written  as: 


m 

dt 


dt  ■ 


dW 


PdQ 


dQ 

-pdq 


dQ 


dW 


dq  ■ 


dq 

{Pi  -  Pi")  dt 


(59) 
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from  which  it  follows  that 
p__^  _  dW 

dQ’  ^  dq' 
dW 

'Hiq,  p,t)+  —  ^  'H\Q,  P,  t). 
ot 


and 

.  _  dm  _  dm  dq  dm  dp 

dQ  dq  dQ  dp  dP 

(60) 

into  the  expression  for  the  velocity: 


(65) 


Inserting  the  second  equation  into  the  third  and  assuming  that 
'H*{Q,  P,  f)  is  simply  a  constant,  we  get  the  famous  Jacobi 
equation: 


■K 


dW 
dq  ’ 


(61) 


whose  solution  furnishes  the  transformation-generating  func¬ 
tion  W.  The  elegant  power  of  the  method  becomes  most  vis¬ 
ible  if  the  constant  'H*  is  set  to  zero.  Under  this  assumption 
the  reduced  two-body  problem  is  easily  resolved.  Starting  with 
the  three  spherical  coordinates  and  their  canonical  momenta  as 
(q,  p),  one  arrives  to  canonically  conjugate  constants  (Q,  P) 
that  coincide  with  the  Delaunay  elements  (27): 

(Qu  Pi)  =  (L,-Mo);  (Q2,  Pi)  =  (G,  -co); 

(G3,  P3)  =  (H,  -O). 

Extension  of  this  approach  to  the  Al-particle  problem  begins 
with  consideration  of  a  disturbed  two-body  setting.  The  num¬ 
ber  of  degrees  of  freedom  is  still  the  same  (three  coordinates  q 
and  three  conjugate  momenta  p),  but  the  initial  Hamiltonian  is 
perturbed: 


d(^  +  m) 

dp 


diPi  +  m) 
- ^ 


(62) 


While  in  (59)-(61)  one  begins  with  the  initial  Hamiltonian  p{ 
and  ends  up  with  TT*  =  0,  the  method  may  be  extended 
to  the  perturbed  setting  by  accepting  that  now  we  start  with 
a  disturbed  initial  Hamiltonian  PI  +  ATY  and  arrive,  through 
the  same  canonical  transformation,  to  an  equally  disturbed 
eventual  Hamiltonian  7T*  -H  ATT  =  A7T.  Plugging  these  new 
Hamiltonians  into  (59)  leads  to  cancellation  of  the  distur¬ 
bance  im  on  the  right-hand  side,  whereafter  one  arrives  to  the 
same  equation  for  W{q,  Q,  f)  as  in  the  unperturbed  case.  Now, 
however,  the  new  canonical  variables  are  no  longer  conserved 
but  obey  the  dynamical  equations: 


dm  ^  _  dm 

dP  ’  “ 


(63) 


Because  the  same  generating  function  is  used  in  the  perturbed 
and  unperturbed  cases,  the  new,  perturbed,  solution  {q,  p)  is  ex¬ 
pressed  through  the  perturbed  “constants”  Q{f)  and  P{f)  in  the 
same  manner  as  the  old,  undisturbed,  q  and  p  were  expressed 
through  the  old  constants  Q  and  P.  This  form-invariance  pro¬ 
vides  the  key  to  the  N-particle  problem:  one  should  choose 
the  transformation-generating  function  W  to  be  additive  over 
the  particles  and  repeat  this  procedure  for  each  of  the  bodies, 
separately. 

Armed  with  this  preparation,  we  can  proceed  to  uncover 
the  implicit  gauge  choice  made  in  using  the  Hamilton- Jacobi 
method  to  derive  evolution  equations  for  the  orbital  elements. 
To  do  this  we  substitute  the  equalities: 


dm  _  dm  dq  dm  dp 

dP  dq  dP  dp  dP 


(64) 


^  dt  dQ^  dP 
This  leads  to: 


(66) 


dq  I  dq  dq  dq  dq\  d  A7T 

'  dq  dp  dq  dp\  d  Ap{ 
[dQ  dP  dP  dQi  dp 


(dm\ 

dp  1 


9  = 


q,  t 


dq 
dt  ’ 


(67) 


where  we  have  taken  into  account  that  the  Jacobian  of  the 
canonical  transformation  is  unity: 


dq  dp  dq  dp  ^ 


(68) 


To  establish  the  link  between  the  regular  variation-of-constsnts 
method  and  the  canonical  treatment,  compare  (67)  with  (41). 
We  see  that  the  symplectic  description  necessarily  imposes  a 
particular  gauge  0  =  5  API  j dp. 

It  can  be  easily  demonstrated  that  this  special  gauge  coin¬ 
cides  with  the  generalised  Lagrange  gauge  (30)  discussed  in 
Sect.  2.2.  To  that  end  one  has  to  compare  the  Hamilton  equa¬ 
tion  for  the  perturbed  Hamiltonian  (19), 


d{P(  +  AP() 

- 


=  P  + 


dAP( 

dp 


(69) 


with  the  definition  of  momentum  from  the  Lagrangian  (17), 


d  (X.(q,  q,  t)  +  AX.(q,  q,  t))  .  d  AX. 

dq  ^  dq 


(70) 


Equating  the  above  two  expressions  immediately  yields: 


0  = 


dAP{\ 

dp  L.t 


(71) 


which  coincides  with  (30).  Thus,  the  transformation  generated 
by  W(q,  Q,  f)  is  canonical  only  if  the  physical  velocity  q  is  split 
in  a  fashion  prescribed  by  (67),  i.e.,  if  (71)  is  fulfilled.  This  is 
exactly  what  our  Theorem  from  Sect.  2.2  says. 

To  summarise,  the  generalised  Lagrange  constraint,  <1)  = 
-d  AXIdq,  is  tacitly  instilled  into  the  Hamilton- Jacobi  method. 
Simply  by  employing  this  method  (at  least,  in  its  straightfor¬ 
ward  form),  we  automatically  fix  the  gauge^.  By  sticking  to  the 
Hamiltonian  description  we  sacrifice  gauge  freedom. 

^  An  explanation  of  this  phenomenon  from  a  different  viewpoint  is 
offered  in  Sect.  6  of  Efroimsky  (2003),  where  the  Delaunay  equations 
are  derived  also  through  a  direct  change  of  variables.  It  turns  out  that 
the  outcome  retains  the  symplectic  form  only  if  an  extra  constraint  is 
imposed  by  hand. 
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Above,  in  Sect.  2.3,  we  established  that  in  the  generalised 
Lagrange  gauge  the  momentum  coincides  with  g.  We  now  can 
get  to  the  same  conclusion  from  (67),  (70)  and  (71): 

d  {£{q,  q,  f)  +  h£{q,  q,  t))  . 

P  =  - -  =^-0  =  0.  (72) 

oq 

Thus,  implementation  of  the  Hamilton- Jacobi  theory  in  ce¬ 
lestial  mechanics  demands  the  orbital  elements  to  osculate  in 
phase  space.  Naturally,  this  demand  reduces  to  that  of  regular 
osculation  in  the  simple  case  of  velocity-independent  AX. 

5.  Conclusions 

In  the  article  thus  far  we  have  studied  the  topic  recently 
raised  in  the  literature:  the  planetary  equations’  internal  sym¬ 
metry  that  stems  from  the  freedom  of  supplementary  condi¬ 
tion’s  choice.  The  necessity  of  making  such  a  choice  con¬ 
strains  the  trajectory  to  a  9-dimensional  submanifold  of 
the  12-dimensional  space  spanned  by  the  orbital  elements  and 
their  time  derivatives.  Similarly  to  the  field  theory,  the  choice 
of  the  constraint  (=the  choice  of  gauge)  is  vastly  ambiguous 
and  reveals  a  hidden  symmetry  instilled  in  the  description  of 
the  A-body  problem  in  the  language  of  orbital  elements. 

We  addressed  the  issue  of  internal  freedom  in  a  sufficiently 
general  setting  where  a  perturbation  to  the  two-body  problem 
depends  not  only  upon  positions  but  also  upon  velocities.  Such 
situations  emerge  when  relativistic  corrections  to  Newton’s 
law  are  taken  into  account  or  when  the  variation-of-constants 
method  is  employed  in  rotating  systems  of  reference. 

We  derived  the  most  general  form  of  the  gauge-invariant 
perturbation  equation  of  celestial  mechanics,  written  in  terms 
of  a  disturbing  force.  Then  we  transformed  it  into  the  most  gen¬ 
eral  gauge-invariant  perturbation  equation  expressed  through 
the  Lagrangian  disturbance. 

Just  as  a  choice  of  an  appropriate  gauge  simplifies  solution 
of  the  equations  of  motion  in  electrodynamics,  an  alternative 
(to  that  of  Lagrange)  choice  of  gauge  in  the  celestial  mechan¬ 
ics  can  simplify  orbit  calculations.  We  provided  one  such  ex¬ 
ample,  a  satellite  orbiting  a  precessing  planet.  In  this  exam¬ 
ple,  the  choice  of  the  generalised  Lagrange  gauge  considerably 
simplifies  matters.  To  achieve  this  simplification,  we  not  only 
exercised  our  right  to  choose  a  convenient  gauge,  but  we  also 
chose  a  preferred  coordinate  system  in  which  to  implement  the 
variation-of-constants  method.  This  interplay  of  the  two  types 
of  freedom  enabled  us  to  eliminate  some  of  the  mathematical 
complications  associated  with  the  inertial  forces.  Not  surpris¬ 
ingly,  it  has  turned  out  to  be  convenient  to  define  the  orbital 
elements  in  the  precessing  frame  of  the  planet;  however,  for 
the  sake  of  mathematical  simplification,  it  also  turned  out  to  be 
beneficial  to  make  these  elements  osculate  in  a  different,  iner¬ 
tial  frame  of  reference. 

We  have  explained  where  the  Lagrange  constraint  tacitly 
enters  the  Hamilton- Jacobi  derivation  of  the  Delaunay  equa¬ 
tions.  This  constraint  turns  out  to  be  an  inseparable  (though  not 
easily  visible)  part  of  the  method:  in  the  case  of  momentum- 
independent  disturbances,  the  A-body  generalisation  of  the 
two-body  Hamilton- Jacobi  technique  is  legitimate  only  if  we 
use  orbital  elements  that  are  osculating.  In  the  situation  where 


the  disturbance  depends  not  only  upon  positions  but  also  upon 
velocities,  another  constraint  (which  we  call  the  “generalised 
Lagrange  constraint”)  turns  out  to  be  stiffiy  embedded  in  the 
Hamilton- Jacobi  development  of  the  problem. 

Unless  a  specific  constraint  (gauge)  is  imposed  by  hand, 
the  planetary  equations  assume  their  general,  gauge-invariant, 
form.  In  the  case  of  a  velocity-independent  disturbance,  any 
gauge  different  from  that  of  Lagrange  drives  the  Delaunay  sys¬ 
tem  away  from  its  symplectic  form.  If  we  permit  the  disturbing 
force  to  depend  also  upon  velocities,  the  Delaunay  equations 
retain  their  canonicity  only  in  the  generalised  Lagrange  gauge. 
Interestingly,  in  this  special  gauge  the  instantaneous  ellipses 
(hyperbolae)  osculate  in  phase  space. 

Briefly  speaking,  A-body  celestial  mechanics,  expressed  in 
terms  of  orbital  elements,  is  a  gauge  theory,  but  it  is  not  strictly 
canonical.  It  becomes  canonical  in  the  generalised  Lagrange 
gauge. 
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Appendix 


In  this  Appendix  we  set  up  an  apparatus  from  which  one  may 
evaluate  the  partial  derivatives  with  respect  to  the  orbital  ele¬ 
ments  of  inertial  terms  that  appear  in  the  planetary  equations 
derived  in  Sects.  3.3  and  3.4.  We  then  show  that  some  of  these 
derivatives  vanish.  Following  that,  we  derive  explicit  expres¬ 
sions  for  each  derivative  of  /r  ■  (/  x  g),  which  provides  a  com¬ 
plete  analytic  evaluation  of  the  rotational  input  in  the  gener¬ 
alised  Lagrange  gauge.  The  topic  is  further  developed  (and  the 
appropriate  generalised  Lagrange  system  of  equations  is  pre¬ 
sented)  in  Efroimsky  &  Goldreich  (2003). 

To  find  the  explicit  form  of  the  dependence  /  =  /(C,,  f),  it 
is  conventional  to  introduce  an  auxiliary  set  of  Cartesian  coor¬ 
dinates  g,  with  an  origin  at  the  gravitating  centre,  and  with  the 
first  two  axes  located  in  the  plane  of  orbit.  The  g  coordinates 
are  easy  to  express  through  the  semimajor  axis  a,  the  eccentric¬ 
ity  e  and  the  eccentric  anomaly  E: 

q\  =  a{co?,E  -  e) ,  q2  =  a^l  -  sinE,  q^  -  0,  (73) 

where  E  itself  is  a  function  of  the  semimajor  axis  a,  the  eccen¬ 
tricity  e,  the  mean  anomaly  at  epoch,  Mq,  and  the  time,  t.  The 
time  dependence  is  realised  through  the  Kepler  equation 

E  -  e  sinE  —  M,  (74) 


where 

M  =  Mo  +  f  At. 

Jto 

The  inertial-frame-related  position  of  the  body  reads: 
r  -  f  (Q,  i,  CO,  a,  e,  Mq;  t)  - 


(75) 


R{Q.,  i,  Lo)  q  (fl,  e,  E(a,  e,  Mq,  t)) , 


(76) 
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R(Q.,  i,  co)  being  the  matrix  of  rotation  from  the  orbital-plane- 
related  coordinate  system  q  to  the  fiducial  frame  (x,  y,  z)  in 
which  the  vector  r  is  defined.  This  rotation  is  parametrised  by 
the  three  Euler  angles:  inclination,  /;  the  longitude  of  the  node, 
Q;  and  the  argument  of  the  pericentre,  a>. 

In  the  unperturbed  two-body  setting  the  velocity  is 
expressed  by 


g  =  —  f(n,  i,  o),  a,  e,  Mq\  t)  = 
ot 


dE 

dt 


a,  e.  Mo 


I  dq 


(77) 


One  can  similarly  calculate  partial  derivatives  of  /  with  respect 
to  Mq\ 

——  /(fi,  i,  u,  a,  e.  Mo;  t)  = 
dMo 


dE  \ 


a,  e,  t 


'  dq 


(78) 


whence  it  becomes  evident  that  df  j dMo  is  parallel  to  g  and, 
hence. 


df  \ 


=  0. 


(79) 


By  a  similar  trick  it  is  possible  to  demonstrate  that  d{f  xg)l dMo 
is  proportional  to  d(f  x  g)ldE  and,  therefore,  to  d(f  x  g)ldt. 
Hence,  this  derivative  vanishes  (because  in  the  two-particle 
case  the  cross  product  /  x  gr  is  an  integral  of  motion).  This 
vanishing  of  d{f  x  g)ldMo,  along  with  (79),  implies: 


fx 


dg  \ 

dMoj^  j^  ^ 


=  0. 


(80) 


In  the  situation  when  the  parameters  C,  are  implemented  by  the 
Delaunay  elements,  a  similar  sequence  of  calculations  leads  to 


gx 


5Mo/q_^  L,  G,H,t 


=  0 


and,  appropriately,  to: 


(81) 


fx 


dg  \ 

dMo]^^  g,h,» 


=  0. 


(82) 


We  can  proceed  much  farther  in  the  generalised  Lagrange 
gauge,  at  least  in  so  far  as  derivatives  of  the  rotational  input 
J  In  —  f  X  g  are  concerned.  (We  remind  that  here  and  ev¬ 
erywhere  n  stands  for  the  reduced  mass,  while  p  denotes  the 
precession  rate.) 

As  we  proved  above,  this  cross  product  is  independent 
of  Mo  and,  hence. 


E  ■ 


d(f  X  g) 

dMo 


(83) 


Since  J  is  orthogonal  to  the  orbit  plane,  it  is  invariant  under 
rotations  of  the  orbit  within  its  plane,  whence 


E  ■ 


d(f  X  g) 
do) 


To  continue,  we  note  that  in  the  two-body  setting  the  ratio  J  fx, 
is  known  to  be  equal  to  ^ja{\  -  e^)  w  where  w  is  a.  unit  vec¬ 
tor  perpendicular  to  the  unperturbed  orbit’s  plane.  Moreover,  in 
planet-associated  noninertial  axes  (x,  y,  z)  with  corresponding 
unit  vectors  (x,  g,  z),  the  normal  to  the  orbit  is  expressed  by: 

w  —  X  sin  i  sin  Q  -  ^  sin  i  cos  -H  z  cos  i.  (85) 


Hence, 

d(fxg)  ^  d(^Ja(l  -  e^)) 

E  ■  - 7 - ■  w - - - 

da  da 

and 


1 

2 


d  if  xg) 


d  ( •\/fl(l  -  e^)) 


de 


Ex  (87) 


where  jx^.  -  Ex  sin  i  sin  Q.  -  jXy  sin  i  cos  Q.  +  jXz  cos  i  is 
the  orthogonal-to-orbit  component  of  the  precession  rate.  The 
remaining  two  derivatives  look: 


and 

d(fxg)  /— t;  dw 

J  ■  — ^  =  V"  0  - ‘■9  *<  ■  -57  = 

f  (1  -  e^)  ^fXx  cos  i  sin  Q.-  jXy  cos  i  cos  O.-  jx^  sin  /j .  (89) 

As  for  the  derivatives  of  a  ■  /,  they  may  be  calculated  di¬ 
rectly  from  the  expression  for  /(fl,  a>,  i,  a,  e.  Mo',  f)  presented 
above.  However,  the  resulting  expressions  are  cumbersome  so 
we  do  not  present  them  here. 
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